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Open quantum system dynamics of random unitary type may in principle be fully undone. Closely 
following the scheme of environment-assisted error correction proposed by Gregoratti and Werner 
[M. Gregoratti and R. F. Werner, J. Mod. Opt. 50 (6), 915-933 (2003)], we explicitly carry out all 
steps needed to invert a phase-damping error on a single qubit. Furthermore, we extend the scheme 
to a mixed-state environment. Surprisingly, we find cases for which the uncorrected state is closer 
to the desired state than any of the corrected ones. 
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I. INTRODUCTION 

Coherent superpositions of quantum states lay the 
foundation for genuinely non-classical phenomena such as 
entanglement or interference of massive particles. With 
growing system size these superpositions become however 
increasingly fragile, thereby impeding future realizations 
of quantum technology. The interaction of the system of 
interest with its surroundings leads to a loss of the inter- 
ference potential — a process usually denoted with deco- 
herence Phase damping (a.k.a. dephasing) denotes 
the case of pure decoherence. Here, the populations re- 
main unaffected; rather, only the coherences itself (i.e., 
the off-diagonal elements of the density matrix) are sub- 
ject to decay. This discrimination of decoherence against 
pure decoherence may for example be adequate in situ- 
ations where the characteristic time describing polariza- 
tion decay is much smaller than the typical dissipation 
time [2]. 

Arguably, one may identify two main mechanisms re- 
sponsible for decoherence. The first and most prominent 
approach rests on the language of open quantum systems 
[3, 3] • Here, the system of interest is seen as being part 
of a larger (often infinite) quantum system, where also 
the environmental degrees of freedom are incorporated. 
The decoherence is then a direct consequence of grow- 
ing correlations between system and environment. It is 
worth noting that these correlations are not necessar- 
ily of quantum nature: the system may in fact decohere 
completely without being entangled to its surroundings 
at all Hi. As a second relevant source of decoherence 
one may consider stochastic fluctuations of ambient fields 
(or "random external fields" 0,I1|)- In fact, fluctuating 
fields were successfully identified as the main source of 
decoherence in ion trap quantum computers, where fiuc- 
tuations are present both in the trapping field and in the 
lasers addressing the individual ions 0. This ensemble- 
based approach may be described in terms of a stochas- 
tic Schrodinger equation, where time evolution is unitary, 
yet stochastic. Usually, dynamics of this type is denoted 
random unitary (RU). 

In this context it is sometimes argued that the latter 
mechanism should be merely seen as a fake decoherence 
process. The argument is partly based on the idea that. 



due to the unitary character of time evolution for individ- 
ual members of the ensemble, the dynamics is in principle 
reversible ^ J,] (the spin echo technique, for example, is 
based on this idea). However, it is known that genuine 
open quantum system dynamics may be reversible, too, 
provided the reduced dynamics is RU and thus indistin- 
guishable from the dynamics caused by fluctuating fields 
|10| . Here, the actual correction procedure is conditional 
on classical information obtained from a measurement 
performed on the system's quantum environment. It thus 
represents an instance of environment- assisted error cor- 
rection. In principle, any phase damping of a single qubit 
or qutrit may be corrected this way, since the dynamics 
is always of RU type [H, V2] . 

The present article provides a thorough assessment of 
all the steps involved in a complete reversal of a given 
phase-damping dynamics on a single qubit. In Sec. 
we discuss the basic tools needed for the environment- 
assisted error correction of RU dynamics. The scheme 
is then applied in Sec. |in]to the example of single-qubit 
phase damping. In addition to the original scheme, where 
the environment is assumed to be in a pure state initially, 
we study a possible extension of the correction scheme in 
case of a mixed-state environment (Sec. IIV[) . 



II. ENVIRONMENT-ASSISTED ERROR 
CORRECTION 

In the theory of quantum error correction one usually 
assumes that there is a certain (low) probability for an 
error acting locally on qubits or gates. Furthermore it 
is assumed that one needs to account for different kinds 
of errors represented by a certain set of error operations. 
The probabilities for different errors to occur are indepen- 
dent of each other and given a priori [Tsj . The correction 
mechanisms are designed such that they can detect and 
correct all errors from the set of possible error operations. 
This requires an encoding of the logical qubits used for 
the actual computation into a larger number of physical 
qubits as can be nicely seen in the Shor code picture [lij . 
In addition, the realization of a unitary gate in the algo- 
rithm requires a translation into a set of unitaries acting 
on the physical qubits. The whole process is described 
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by the beautiful theory of fault tolerant quantum compu- 
tation Unfortunately as pointed out in, e.g.fis!], the 
scaling of the number of physical qubits is, though only 
polylogarithmic in the asymptotic limit, too severe even 
for a very small number of computational qubits to per- 
mit an experimental realization by current means. Some- 
times, being able to perform measurements on the envi- 
ronment of the qubits allows to recover quantum states. 
As an example, in [l5| a continuous monitoring of the de- 
cay processes provides the information for a conditional 
recovery operation. 

Environment-assisted error correction as proposed in 
is a different approach to quantum error correction. 
Here, the correction is based on classical information ob- 
tained from a measurement on the system's environment. 
The scheme allows for a complete correction of RU dy- 
namics, provided the (quantum) environment is initially 
in a pure state. We closely follow the ideas presented in 
[lOj : we are, however, interested in carrying out the pro- 
cedure lined out in the aforementioned work for a feasible 
realization of an open quantum system. In the following 
we will present a scheme that allows the explicit construc- 
tion of an observable for such a measurement. Our results 
will provide a correction scheme for a phase-damping in- 
teraction between a qubit and an initially pure environ- 
ment of finite dimension. 



where Pa > 0,J2aPa = 1- The Kraus operators are thus 
unitary up to normalization, Ka — ^JVaUa- Obviously, 
a RU channel is unital, that is, leaving the completely 
mixed state invariant, <J>ru(1) = 1- 

Phase damping (or dephasing) is the case of pure 
decoherence, where, in a fixed basis {In)} (the phase- 
damping basis), no population transfer occurs. Accord- 
ing to this basis of "preferred states" , the projectors are 
constants of the motion. In such a case the Hamilto- 
nian describing the system-environment coupling may be 
diagonalized with respect to the phase-damping basis: 
H = 'Ylin l^-X"-! ® Here, the relative Hamiltonians 

h„ act on the environment, only. Assuming the environ- 
ment to start in a pure state, IV'o): the phase-damping 
dynamics is fully described in terms of the overlap of the 
relative states IV'n) := e~*"*|'(/'o), that is 

Pmn = (V'n|l/'m) Pm«- (5) 

Leaving the diagonal elements intact, phase-damping 
channels are unital by definition. In their article 
Landau and Streater show that a phase-damping chan- 
nel acting on a single qubit or qutrit may always be de- 
composed into a RU decomposition, Eq. In principle, 
phase-damping errors on systems of dimension 2 or 3 may 
thus be completely undone. 



A. Quantum channels 

In the context of open quantum systems, a quantum 
channel describes the time-evolution of a quantum sys- 
tem arising from the joint unitary evolution of system 
plus environment. Under the usual assumption of van- 
ishing initial correlations, system and environment are 
initially described by a product state p® pE- The chan- 
nel is then obtained by tracing out the environmental 
degrees of freedom 

$(p)=tr£(e-'"V<»P£e'H*). (1) 

The total Hamiltonian H includes the interaction be- 
tween system and environment. According to \v?\ . 
every such channel has a nonunique decomposition 

p' Hp) = K^PKt, = E L/3PLt , (2) 

where two corresponding sets of so-called Kraus opera- 
tors Ka,Lp are related via a unitary matrix V = (wa^) 
pj , so that 

KQ=^Ua^L/3. (3) 

13 

In case of RU dynamics it is possible to find a decom- 
position of the channel into unitary operators, 

"^lixjip) ^^Pc^^cP^l, (4) 



B. The correction scheme 

If the initial state of the environment is pure, pe = 
l''/'o)('0o|j a certain Kraus decomposition is selected by 
choosing a basis {|x/3)}/3 of "He- Upon insertion into 
Eq. H]), this leads to 

Hp) - Y.^xp\e-''''\i^,)p{M^'''\x0) 

where we identify L/j = (X/sl^^'^'I'^o)- 

In principle, this identification allows to select the dy- 
namics where only a single term of the operator sum ([2]) 
applies: 

L^pL^ - ivE {e-'^'p ® PE(^™{t Pp)) , (7) 

with the projector = |x^)(x,3|. The unitary equiv- 
alence between different decompositions allows to write 

Ka = (Ma|e~™|l/'o) with \pa) = E/3 Va/3lX/3)- 

The principle of correction of a RU channel ^rjj is 
now straightforward: it relies on the identification of an 
appropriate basis \pa) corresponding to a RU decompo- 
sition of the channel. The projectors Qa — \pa){Pa\ may 
then be used to single out the subnormalized sub-state 

Pa = PaUapUj, = tlE (e~'"*p ® PBe'H*(l (g, Qa)) , 

which, upon normalization, is unitarily related to the ini- 
tial state p. 
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For a suitable observable O on He with projectors 
on non-degenerate eigenspaces, O — ^aQa with 
tr(QQQ^) = Sa^Qa, a measurement outcome of Aq, al- 
lows to unambiguously discriminate between sub-states 
p'^. The perfect recovery of the initial state is achieved 
by UJ,/9qUq. Gregoratti and Werner show in [l3] that 
RU channels are the only ones thus allowing for a perfect 
correction. 



III. EXPLICIT IMPLEMENTATION: 
CORRECTION OF SINGLE QUBIT PHASE 
DAMPING 

In practice, the correction scheme faces several imped- 
iments. First of all, no simple method is known how to 
decide whether a given channel has a RU decomposition. 
Second, even if such a decomposition is possible, the sin- 
gle unitaries Uq. have to be known in detail. Once these 
hurdles are overcome, one has to find explicit expressions 
for all Qq,. The realization of this last step is discussed in 
the original formulation [loj . Here, the unitary equiva- 
lence between equivalent Kraus decompositions, Eq. ([3]), 
plays a central role. 

Now we have all ingredients to construct a Hamiltonian 
generating a RU channel on a qubit and to proceed with 
the construction of the observable. For a phase-damping 
channel on a qubit there are only two relative Hamil- 
tonians hi, h.2 leading to the relative states \^i{t)) = 
e""'i*|V'o), |V'2(i)) =e""'2*|V'o) on the environment. The 
action of the phase-damping channel has a very simple 
form in terms of their overlap, C{t) := {ip2{t)\ipi{t)), 



P = 



Pll Cpi2 

Cp2i P22 



(8) 



We have suppressed the time dependence in the equation 
above and will from now on only consider quantities for 
a fixed value of t. We have already seen that using the 
basis {|X;g)};3 oH-Le to explicitly perform the trace allows 
us to identify hp — (x^i|e '^*l'0o) so that we can write 



L« = 



(x/jIV-i) 

(X^IV'2) 



(9) 



A. RU decomposition 

In order to obtain the RU decomposition we com- 
pute the dynamical or Choi matrix of our phase-damping 
channel. For n — dimHs, '■ 'Hs — >■ "Hs we can identify 
Hs ^ C" and the Choi matrix is defined as the ri^ x 
matrix containing the action of $ on all n x n matrices 

E,;i = £ 



that form a basis of A^(n x n, C), 



1 C" 





C 1 




(10) 



Diagonalization of this 4x4 matrix yields the nonzero 
eigenvalues 



7l( 



^,0,0,1)^ 



• Ai = 1 - |C|: VA, - T7^\~jc] 
. A2 = 1 + |C|: -;).( 1^,0, 0,1)^ 



with corresponding eigenvectors v^^ , vx^ ■ According to 
a central result in [18| we can obtain Kraus operators by 
rearranging the eigenvectors of the above Choi matrix 
into 2x2 matrices, resulting in the decomposition 



1-\C\ 



l + \C\ 



c 
'\c\ 


c 




Jc\ " 
1 

Jc] " 
1 



(11) 
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E 



K^pKl 



Note that the Kraus operators of the decomposition in 
(ITU) are trivially unitary up to a scaling factor, KaK^^ = 
K\Ka = PqI, a = 1, 2. It is easily verified that this RU 
decomposition recovers Eq. ([8]). 



B. Finding the correction observable 

Consider now the situation where the actual environ- 
ment is of dimension n. As a consequence the set 
contains n elements and according to Eq. ((S]) we know 
that there exists a unitary matrix V that relates {L^j/j 
to the RU decomposition {Kq}q in pT|) . if the latter is 
extended by n — 2 zero matrices K^ = ^/p^a, o: = 1,2, 
Kq = 0(2) J a — S, . . . ,n. We will now show in detail how 
it is possible to obtain V from the two sets of Kraus oper- 
ators {Kq}q and {L/3}/3. Due to their diagonal character 
the Kraus operators may be rewritten in terms of vectors 
Ka — {Kii, K22)- The unitary equivalence, Eq. ([3|), thus 
translates to the following linear system for the rows Vq 
of V, 



L2 L2 



f Val \ 
V Van J 



(12) 



where in addition all double indices are suppressed. For 
a = 1 , 2 (fT2|) is an inhomogeneous system Av^ = bo, 
with inhomogeneities obtained from Ki, K2 given by 



bi 
b2 



c 
\c\ 

1 



l+lgl I l§l 

2 \ 1 



(13) 
(14) 



For a = 3, . . . , n one has to fulfill the homogeneous sys- 
tem Avq = 0. In addition, unitarity of V requires that 
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(vq|v^) = dap- In other words we need to find n — 2 
orthonormal vectors Vq, S ker(A) and two orthonor- 
mal vectors Vq, e ker(A)-'- satisfying the inhomogeneous 
equations. In the following we will show that the singular 
value decomposition of A will provide all of the required 
solutions Vq,. 

The 2 X n matrix A has singular value decomposi- 
tion A = UEW^, with U a unitary 2x2 matrix with 
columns spanning im(A) , S a 2 x n diagonal matrix con- 
taining the singular values of A, and W a n x n uni- 
tary matrix. If r = rank(A), the first r columns of W 
form an ONB of ker(M)^ and the last n — r columns 
form an ONB of ker(A). Using the completeness re- 
lation X]/3=i \xp){xp\ = 1 it is straightforward to show 

^ '^^ ==U(I]Et)ut. Diagonalizing AAt 

l-\C\ 
1 + |C| 



that AAt = 
results in AA^ 



U 



U\ with eigen- 



vectors Ui = ^ f I'^l 

V2 \ 1 



c 



forming the 



columns of U. The singular values of A can be read off as 
the square root of the eigenvalues of AA'^, Ai — ^1 — |C| 
and A2 = \/l -I- |C|. Now let be the columns of W. 
Using the singular value decomposition of A it is now 
obvious that Awi = AiUi = bi, Aw2 — A2U2 = b2 and 
Awi — for alH — 3, . . . , n. This means that the singular 
value decomposition of A provides the desired solutions 
Vq to the linear system (fT2|) via the columns of the 
unitary matrix W. Thus we can finally conclude that the 
unitary matrix V relating the RU decomposition of 
to the decomposition with respect to the basis {|x)^}/3 is 
given by 

V = W^. 

We can transform the basis {|x),3}^ to the desired mea- 
surement basis {|/iQ)}Q via 




(b) 

Figure 1: (Color Online) Phase-damping and correction 
for an initially pure environment, (a) The pure initial 
state p of the system (black arrow) is mapped to a 
mixed state $t(p) (thick, black line). The 
phase-damping interaction leaves invariant the z-value 
of the Bloch vector (the time-evolved states lie on a 
plane of constant z). (b) The corrected state 
p = Ra{pa) after a measurement of O with outcome Aq. 

and application of the appropriate correcting channel 
Ra ■ The corrected state p is identical to the initial state 
p of (a) for all of the above channels $4. 



(b) displays the corrected state for each of the states 
^t{p) after a measurement of the correction observable 
O (which in general will be different for different channels 
<&() and the application of the correction procedure Ra 
(which depends on the measurement outcome a as well 
as the channel $4). The correction procedure recovers 
the initial state p for each channel <I>t. One can think of 
the combined action of measurement and correction pro- 
cedure as the process that mediates between Fig. [T] (a) 
and (b). 

The phase-damping channel <^t is completely deter- 
mined by the two relative states \'4)i{t)), |V'2(^))- Their 
unitary time evolution is displayed in Fig. [5] (a). The 
states \p)a spanning the eigenspaces of the correction 
observable O do also depend on the time parameter t la- 
beling the different phase damping channels $4. Their 
time-evolution is displayed in Fig.[2](b). 



The projectors Qq for the observable O on He are given 

by Qq = \p,a){Pa\- 

Now the general correction scheme consists of the fol- 
lowing steps: i) Performe a measurement on the environ- 
ment using the observable O = J2a=i ^aQa, with Ai 7^ 
A2 resulting in a post measurement state pa = UapUj, 
for all possible measurement outcomes Aq.. ii) Apply the 
unitary transformation Vl^paVa conditional on the out- 
come of the measurement of O. The probability to get 
Aq is p(Aq) = Pa, a = 1,2 and p(Aq) — otherwise. 

The following sequence of figures intends to provide 
further insight into the mechanism of the correction 
scheme. For simplicity we consider a two dimensional 
environment so that we can represent state vectors of sys- 
tem and environment in the Bloch-sphere picture. Fig. [1] 
(a) shows a sequence of phase-damping channels $4 act- 
ing on the pure initial system state p. The adjacent Fig.[T] 



IV. MIXED ENVIRONMENTS 

Until now it was assumed that the environment is ini- 
tially in the pure state {ipo}- In this section we want 
to investigate the case of a mixed-state environment. 
We apply modified versions of the correction scheme 
for pure environments and explicitly show that none 
of these protocols leads to a satisfying correction of 
the phase-damping channel. For a mixed environment 
there is no unique correspondence between the projec- 
tion Qq on the environment Hilbert space and a single 
Kraus operator Kq as in the case of a pure initial en- 
vironment. We assume that our environment is two- 
dimensional and that its initial state is given by the 
mixture pE — w\tpo){4>o\ + (1 — w)\iPq){iI'o\- For con- 
creteness think of a low but not zero temperature envi- 
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Figure 2: (Color online) (a) Unitary evolution of the 
relative states \i}i{t)) (red) and |V'2(0) (blue). The 
initial state |-i/;o) is indicated by the black arrow lying on 
the intersection of the cones, (b) Correction observables 
El (red) and E2 (blue) for different channels $4. The 
two observables span orthogonal subspaces of He (the 
corresponding Bloch vectors point to opposite points). 



ronment such that w is less but close to 1. The cou- 
pling is realized via the following Hamiltonian, H = 
ka^ ^af + l^Ta, a — J^^^i '^i^i the vector of Pauli-spin 
matrices and F G R^. The Hamiltonian can be written 
as 2 l*2:)('^2l ^ hi, where |s*) are the eigenstates of 
af and hi = kaf + Ta^ , h2 = -kaf + Ta^ . The chan- 
nel $ generated by H according to Eq. ^ preserves the 
diagonal elements of p in the |s^) basis. For pE mixed the 
phase-damping interaction generates two additional rel 



ative states IV'i 



-ihi* 



IVo^) and l^Pi) 



He- The action of $ is similar to the pure case, Eq. ([S|), 
but C is replaced by wC -I- (1 — w)C^ . C"*- is again given 
by an overlap of relative states, C""" = (V'2"IV'i~)- 



A. Correction schemes for mixed environments 

In order to perform a correction we treat the channel as 
if it was generated by the pure initial environment state 
iV'o) (or IV'o"))- For jV'o) we calculate the correspond- 
ing projectors Qi, Q2, the observable O = 2 ^i^i^ 
7^ ^2, for a measurement on the environment, prob- 
abilities p(Ai), p{X2) of the RU decomposition and uni- 
tary operators Ui, U2 for the correction operation. We 
will call the states resulting from this correction pro- 
cedure pa,c- Similarly for we calculate Pi, P2, 
O-"- = I]j=i,2 Ml 7^ ri, r2 and Vi, V2 and caU 
the corrected states aa,c (the possible correction path- 
ways for a mixed state are displayed in Fig. [3]). A mea- 
surement of O (O"*") and a successive outcome dependent 
correction (V^) will not perfectly correct $ since we 
have neglected the fact that p_E is a mixture. 

Several correction protocols can be considered. First, 
the corrector might simply ignore the fact that there is 
a small admixture of IV'o')- this case, the state after 
correction is given by p(Ai)pi^c + p{X2)p2,c ='■ Pc- Sec- 
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Figure 3: Schematic picture of the possible correction 
schemes for a mixed environment (see text for details). 



ond, as a matter of keeping an open mind the correc- 
tor might want to mirror the mixed-state probabilities 
for the two states in his choice of correction scheme. In 
other words, he performs the measurement on He and 
the successive correction with a probability of w as if p^; 
was the pure state |'0o) and with a probability of 1 — w 
as if pe was \4>q). We call the state resulting from an 
ensemble average over states corrected in this manner 
Pc = wpc + (1 - w)ac- 



B. The error of the correction 

In order to quantify the quality of the correction we 
measure the distance between the corrected states Pa,c, 
<ya,c, Pc and Pc and the initial state p, and compare it 
to the distance between the uncorrected state $(p) and 



in the initial state. The trace norm 0, 



\P\\ - 2' 

is a suitable distance measure in the set of states that 
will later allow us to derive strict bounds analytically. 
In Fig. 13] we show results of the correction protocol for a 
small admixture of \tpQ) to \ijjo) in Pe {w = 0.9). Clearly, 
applying any of the correction procedures does not nec- 
essarily lead to an improvement. Indeed, we can see 
that even for a small admixture there are phase-damping 
channels $t (corresponding to interaction times t) for 
which the uncorrected state is closer to the desired state 
than any of the corrected ones. 

We observe from Fig. that the correction scheme is 
likely to worsen the state in cases where the influence of 
the channel is small, i.e., ||p— <I>(p)|j — ^ 0. In order to shed 
light on this worst case phenomenon we derive explicit 
expressions for the trace distances. Channels that do not 
significantly alter the input state can be characterized by 
assuming C = 1 — ie, = 1 + ie with O(e^) ~ 0, see 
App.EI For such channels it is possible to find analytical 
expressions for the trace distances. An explicit derivation 
can also be found in App. [B] For pE mixed, w ^ 1 
they are given (up to terms of ©(e^)) by ||p — $(p)|| = 
2|pi2|h-i|e, ||p-pi,,|| =2|pi2||l + ie|, ||p - P2,c|| = 0, 
and Hp — Pell = 2|pi2||l — w||l -f ie\. We can see that 

||p-Pi,c|| > Hp -Pell > ||p-$(p)|| > I|p-P2,c||. 
Furthermore, we find ||p — Pc|| = 2|pi2||u'(l — w) + ^| and. 
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Figure 4: Trace distance as a measure for the quality of 
the corrected state. Different channels $t arising from 
H and an initial state Pe containing a small admixture 

of {w = 0.9) are corrected using the correction 
scheme for pure states (details see text). Clearly, there 
are cases where the uncorrected state ^t{p) (solid line) 
is closer to the desired state than the results from any 
of the correction schemes. 

hence, 

Hp- Pell > ||p-$(p)||. 

This allows to conclude that for channels which do not 
significantly change the initial state p the uncorrected 
state $(p) is a better approximation to the initial state 
than the corrected states pc, Pc- 

The inequality also shows that an improvement can be 
obtained by selecting only the states p2,c corresponding 
to the very rare outcome A2. Note, however, that for the 
corrector it is a priori impossible to judge which correc- 
tion scheme works best. 



How to correct a phase-damping error arising from 
an environment starting in a mixed state is a delicate 
task. First, a sensible correction protocol needs to be 
defined. We propose two rather straightforward schemes 
and study their ability to correct. We observe that in gen- 
eral the correction procedure for pure environments is no 
longer successful. Indeed, there are cases for which the 
uncorrected state is closer to the desired state than any 
of the corrected ones. The surprising conclusion is that 
even small admixtures to the initial state of the environ- 
ment renders the success of the correction undetermined. 
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Appendix A: The approximation for \C'\ ~ 1 

From the above discussion we know that requiring the 
phase-damping channel to be close to the identity is equal 
to saying that wC + (1 — w)C^ w 1. Expanding the 2x2 
Hamiltonian h^ in terms of Pauli matrices h = Y^^=i li'^i 
and using the closed expression for the unitary operator, 
e^'''* — cos art — isinatncr, a — vSTt? ' n = ^ we 
find that Re(C) = Re(C^) and Im(C) = -lm(C-L). This 
means that for such channels Rc(C)-|-i(2u'— l)Im(C) w 1 
and thus C = 1-ie, = 1+ie with 0{e^) 0. Utilizing 
the closed expression for e~''^'=* we can furthermore show 
that there are times for which the above requirements on 
C, are simultaneously fulfilled. 



V. CONCLUSION 

A phase-damping error of a single qubit arising from 
the interaction with a pure and mixed finite-dimensional 
quantum environment is considered. For the pure en- 
vironment, we apply the correction scheme proposed in 
[lOj using knowledge about the processes governing the 
full system dynamics and the initial state of the environ- 
ment. The correction is based on the RU decomposition 
for every such channel from a spectral decomposition of 
the channels Choi matrix. We are able to explicitly con- 
struct the unitary matrix relating a decomposition with 
respect to an arbitrary basis of the environment to the 
RU decomposition. The dynamics of such an open quan- 
tum system can effectively be described using an envi- 
ronment of only two dimensions spanned by the relative 
states li/ji) and \ip2)- 



Appendix B: Analytical expressions for the trace 
distance 



In order to derive analytical expressions for the trace 
distances we express the post measurement states for a 
measurement of O, pa, and of O"*", aa, in terms of the 
following channels. We define RU channels: ^Ru.a{p) = 
tris (Up(8) |i/'o)(V'o|UU(g)Qa) = PaVapVl, "^ruAp) = 
tri; (Up0 |i/'^)(V'^|UU (g)Pa) = r^Yapyi and auxil- 
iary channels ^err,aip) = trs (Up (g) |V'o)(V'o|U'fl (g) Qq), 

so that we 



«'err,a(p) = tv E {V p ^ \^i^){lP^ \W t ® F . 



can express pa 



trE(u)*R[/,Q(p) + (l-to)*err.a(p)) 



and simi- 



larly (Tq, replacing $ with and interchanging w and 1 — 
w. One can now see that the error of the correction pa^c = 
UjjpaUa comes from the term {1 — w)'Ul^^err,a{p)V a and 



7 



that pa,c contains the desired state p as wp^p. A sim- 
ilar statement is also true for (7a, c- the following we 
will only discuss the states pa.c but a similar argument 
applies to aa,c- Using tr{^Ru,a{p) + ^err,a{p)) = 1 we 
can verify that tr ($err,a(p)) = 1 — Pa- Observing that 
p(Ai) = wtr^Ru,aip) + (1 - w)tv^err,a{p) we Can use 
the above result to write down an analytical expression 
for the probabilities of the measurement outcomes 
of O, p(Aa) = wpa + (1 — w){l — Pa). In thc case 
that C = 1 — ie we can approximate the probabilities 
from the RU decomposition as pi ~ 0, P2 ~ 1- Us- 
ing thc fact that for positive matrices the trace norm 
equals the trace, tr(>l) = trVAtA = ||A||nMc and that 
for finite dimensional matrices the trace norm is equiv- 
alent to the maximum norm |lA||oo = max|ajj| we 

can conclude that ^ru,i{p) — and ^err.2{p) = 

©(e^). Furthermore we can express thc unitary matri- 
ces from the RU decomposition as Ui = ^ ^ 



1 



and Uo = 



1-ie 
1 



corrected states as pi c = ul* 

err. ,i{p)Vi and p2,c = 
can already conclude that p2,c = 
p + 0{e^). In order to obtain an expression for $err,i(p) 

we use that ^{p) = Y,a (.W^RuAp) + (1 - w)^'err,a(p)) 

and ||4'err,2(p)|| = 1 — P2 = SO that we can ap- 
proximate = (1 — w)$err,i(/9) + w'^RUfiip) and fi- 
nally express $err,i in terms of <I> and The fi- 
nal expression for the corrected state pi^c is then given 
pii -(1 -I- 2ie)pi2 
-(l-2ie)p2i P22 
straightforward to calculate \\p — pi,c\\ ~ 2|pi2||l -l- ie], 
=2|pi2||w- \\e, and \\p - P2A\ =0- For the 
corrected state pc we obtain Pc|| = 2|pi2||l— 
Performing a similar calculation for (Ja,c and p{pa) one 
can obtain the following expression for thc distance be- 
tween initial state and the corrected state pc- \\p— Pc\\ = 
2\pi2\\w{l - w) + \\. 



by Pi.c = 



It is now 



This allows to express the 
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